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Abstract Fifteen years after its discovery, the physics of transient seismic emis- 
sion in flares remains largely mysterious. An early hypothesis by its discoverers 
proposed that these "sunquakes" are the signature of a shock driven by "thick- 
target heating" of the flaring chromosphere. H-a observations show evidence 
for such a shock in the low chromosphere during a flare. However, simulations 
of shocks driven by impulsive chromospheric heating show withering radiative 
losses as the shock proceeds downward. The compression of the shocked gas not 
only heats it but increases its density, making it more radiative. So, radiative 
losses increase radically with the strength of the shock. This has introduced 
doubt that sufficient energy from such a shock can penetrate into the solar 
interior to match that indicated by the helioseismic signatures. The purpose 
of this paper is to point out that simulations of acoustic transients driven 
by impulsive heating have yet to account for magnetic fields characteristic of 
transient-seismic-source environments. We show this to be a critically important 
factor that has a major impact on the seismic flux conducted into the solar 
interior. A strong horizontal magnetic field, for example, greatly increases the 
compressional modulus of the chromospheric medium. This greatly reduces com- 
pression of the gas, hence the radiative losses as the transient passes through 
it. This could explain the strong affinity of seismic sources to regions of strong, 
highly inclined penumbra! magnetic fields, including the neutral lines separating 
opposing polarities in (5-configuration sunspots. The basic point, then, is that the 
role of inclined magnetic fields is fundamental to our understanding of the role 
of impulsive heating in transient seismic emission. Obliquely inclined magnetic 
fields will complicate simulations of impulsive heating considerably. However, 
horizontal magnetic fields, as a preliminary control simulation, can be incorpo- 
rated into standard f-D thick-target-heating simulations with a relatively simple 
adaptation of existing HD codes. 
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1. Introduction 



The first known instance of seismic transient emission, also known as a "sun- 
quake," was that discovered by IKosovichev and Zharkoval (11998^ . emanating 
from the X2.6-class flare of 1996 July 9 in the declining phase of Cycle 22. Seven 
years passed, with some major acoustically quiet flares, before Donea and Lindseyj 
(j2005p found seismic emission from multiple sources in the flares of 2003 Octo- 
ber 28-29. This led to the subsequent discovery of about a dozen instances of 
transient seismic emission in the declining phase of Cycle 23, some associated 
with relatively weak flares. Transient seismic emission is now known to be a 
moderately common phenomenon. Nevertheless, something about the mechanics 
behind this phenomenon continues to elude us. 

The significant manifestation of transient seismic emission is a pattern of 
outgoing ripples seen in Doppler-helioseismic observations 15-45 min follow- 
ing the impulsive phase of a flare. Computational seismic holography applied 
by |Donea, Braun, and Lindsey] (|1999p . |Donea and Lindsey] (I2005P and others, 
showed compact sources up to ~12 Mm across at sites of hard X-ray emission, 
and sometimes 7-ray emission. Estimates of the seismic energy, Es, emitted by 
a single compact source, based upon signatures derived from seismic holography 
( [Donea, Braun, and Lindsey, 1999||Donea and Lindsey, 2005^ , are frequently in 
the range 10^''-10^^ erg0 The vertical momentum carried by a seismic transient 
at the solar surface is estimated at Es/c, where c is the sound speed. Taking c 
to be '^lO km s^^ renders seismic momenta in the range 10^^-10^'^ dyne s. 

In large flares, there may be more than one such source, in which case the total 
seismic energy emitted can be correspondingly greater. This acoustic energy is 
only a small fraction of the total energy ( Emslie et al., 2012 ) estimated for a 
flare, generally a fraction of a percent, and for large flares often less than 10^"* 
of the total. 

Attempts to understand transient seismic emission have focused upon three 
mechanisms: 



1.1. Impulsive Heating of the Chromosphere 

IKosovichev and Zharkoval (|1998p attributed transient seismic emission to a pres- 
sure wave driven by "thick-target heating" of the chromosphere by energetic 
electrons. A good deal of what we understood of this mechanism at the time 
came out of 1-D numerical HD simulations (Fisher, Canfield and McClymont 
1985a, b,c, thenceforth FCM) applied to pre-flare models of plage chromospheres 



^Computational seismic holograpliy is based upon the extrapolation of the acoustic field 
from a surface region several thousand km from a source region back to the source region 
itself. Estimates of the energy flux from this signature are the acoustic analogy of pho- 
tometry in the electromagnetic spectrum. For an example of this application, we refer to 
lAlvarado Gomez et al.l II2012II . 
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and coronae. Energetic electrons accelerated in the corona heat the upper chro- 
mosphere to something like coronal temperatures. The initial pressure of the 
medium, 10-100 dyne cm~^ is increased by about two orders of magnitude. 
The resulting increase in pressure drives a shock downward into the medium 
beneath the heated layer, carrying a momentum of the order of lO^'^ erg cm 
(jKosovichev and Zharkoval I1995p . This is essentially balanced by the upward 
momentum of the heated layer, which explodes upward into the pre-flare corona. 
IZarro and Canfieldl (I1989P found this shock model consistent with H-a emission 
spectra typical of M-class flares. 

More recent consideration to the role of impulsive chromospheric heating in 
transient seismic emission is given bv lZharkoval (|2008p and lZharkov et aL](|201ip . 

The hypothesis that transient seismic emission is driven by impulsive chromo- 
spheric heating has the problem that severe radiative losses deplete the transient 
as it passes through the underlying low chromosphere. These concerns were 
first communicated to us by Fisher (2006, private communication). They are 
mentioned by Hudson, Fisher and Welsch] (|2008p and treated in further detail 
bv IFisher et a .1 (|2012p . This conclusion is reinforced by FCM-style simulations 
bv lAllred et al.l (|2005p . and, with improved treatment of radiative transfer, by 
Allred (private communication, 2012). Even in acoustically active flares the 
signature of heating is often evident in regions from which no seismic emission 
appears to emanate. This includes acoustically inactive flares, from which no 
transient seismic emission at all is evident. 

Our understanding of the dynamics and geometry of atmospheric heating has 
changed dramatically since FCM. The magnetic structure of penumbrae, from 
which transient seismic emission generally appears to emanate, is highly fila- 
mentary, as opposed to the plane-parallel models of FCM. IFletcher et al.l (I2008P 
suggest that the high-energy electrons that serve as the primary agent for the evi- 
dent heating could be accelerated much deeper than in the models of FCM, i.e., in 
the chromosphere instead of the corona. Indeed, recent 5'£'0/STERE0/RHSSI 
observations by IMartinez et al.l (I2012p using parallax to determine the source 
heights of white light (WL) and hard X-ray (HXR) emission associated with 
atmospheric heating indicate that this heating is significantly deeper than in the 
models of FCM. 

1.2. Backwarming of the Photosphere 

Based upon some of these developments, |Donea, Braun, and Lindsey] (I1999P 
and IMoradi et al.l (|2007p proposed that heating of the photosphere by visible 
and near-UV continuum radiation from the heated chromosphere could drive a 
seismic transient (see also [Lindsey and Donea[ |2008p . In this model, a pressure 
transient is launched at the base of the photosphere, which most efficiently ab- 
sorbs the continuum radiation from the chromosphere, and proceeds downward 
(and upward) from thence. Once the downward component of the transient 
penetrates into the solar interior, further radiative losses are blocked by the 
high opacity of ionized hydrogen. This mechanism, known as "backwarming," 
was suggested by sources of seismic transient emission that coincided closely 
with transient white-light flare kernels ( [Lindsey and Donea[ 12005)) . 
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Instances have since been discovered of strong seismic emission from regions 
that show no appreciable excess in visible continuum radiation integrated over 
the neighborhood of the seismic source. An example is the X2.2-class flare 
on 2011 February 15. Integrated over a circular disk whose area is 90 Mm^, 
encompassing the seismic source, these fluctuations are only 0.12% of the quiet- 
Sun intensity. For comparison, a similar average taken over the seismic source 
of the white-light flare of 2003 October 29 found a continuum emission 10% 
in excess of the quiet-Sun intensity ( |Donea and Lindsey, 2005 1, i.e., about 80 
times greater than in the flare SOL2011-02-15 ( Alvarado Gomez et al., 2012[ ). 
According to our rough understanding of how back- warming could prospectively 
work as a seismic source ( [Lindsey and Donea, 2008 ) , a compression wave driven 
by heating commensurate with this degree of fluctuation in continuum intensity 
could explain barely 2% of the the seismic energy that appears to emanate from 
the flare SOL2011-02-15. Because of this, it is evident that some other mechanism 
must be in operation to explain the seismic transient released by this flare. 

Variations to the back- warming theme draw upon energetic protons ( Zharkova, 2008| .[ 
Protons with energies of 200 MeV could penetrate to the base of the photosphere, 
heating it directly. Protons of these energies interact with the nuclei of elements 
such as oxygen, iron and magnesium ( [Najita and Orrall, 1970"] ) to excite 7 rays. 
Such 7 rays are indeed observed in some seismically active flares. However, we 
now know of instances (IDonea et al.[ I2006P of strong seismic emission from flares 
from which no 7 radiation was detected. 



1.3. Lorentz-Force Transients 



Based largely upon the foregoing assessments, [Hudson, Fisher and Welsch| ([^0051 1 
and IFisher et al.l (|2012p proposed that transient seismic emission is the result of 
Lorentz-force transients, caused by the sudden re-configuration of the magnetic 
field in the impulsive phase of the flare. 

Magnetic signatures roughly consistent with Lorentz-force transients driving 
transient seismic emission consistent with helioscismic signatures are reported 
by IDonea et al.l ISudol et al.l [Wang et al.[ ([^00^ [Wang and Uu\ 

PUTU)) . IPetrie and Sudoll ([^giUj) . iPetnel f^UT^ and others. These signatures ap- 
pear to be encumbered by two significant qualifications: 



• Magnetic signatures representing highly perturbed radiative environments 
such as are typical of flares are notoriously unreliable, often showing mo- 
mentary reversals in the field direction inconsistent with expectations based 
upon MHD in a highly conducting medium. 

• Such magnetic signatures are common in acoustically inactive flares, as 
they are in acoustically active flares in wide-spread regions from which no 
detectable seismic emission emanates. 



lAlvarado Gomez et al.l(|2012p studied the seismically active flare SOL2011-02- 
15, concentrating on the magnetic signature in the source region of the strongest 
of two sources (Figure 1). Based upon line-of-sight magnetic signatures immedi- 
ately before the impulsive phase and 10-15 minutes after, they offer a somewhat 
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discouraging appraisal of the potential of the local Lorentz force transient as 
a driver for the seismic emission observed. However, they do not rule it out. 
Further consideration to the role of Lorentz-force transients in this flare is given 
bv lZharkov et aL]p011|) . based largely upon the apparent magnetic connectivity 
between the primary acoustic source and a second one, which they themselves 
discovered. 

AIA observations show strong 160 nm emission from the end of a ribbon 
^2.8 Mm SSE of the seismic source in the early impulsive phase (middle- 
right panel of Figure 1). Within seconds, emission in this region is obscured 
by saturation/overflow from strong, extended emission at large. Hinode/SOT 
observations in the Call K-linc show a very similar ribbon that sweeps northward 
through the acoustic source region in the minute or so following the time of 
the 160 nm snapshot. This suggests strong heating of the chromosphere, that 
might drive a seismic transient — if the resulting shock could somehow escape 
the radiative losses described in Hl.ll above. lAlvarado Gomez et al.l (|2012p re- 
port HXR emission whose RHESSI signature encompasses the seismic source 
region, which would connect the heat source to energetic electrons!! Like the 
magnetic signatures, both the 160-nm and the HXR sources are stronger in 
other far-removed regions, from which no seismic emission is detected. 

2. Regional Selectivity 

In summary, then, the signatures indicating impulsive chromospheric heating 
and Lorentz-force transients appear to share one conspicuous quality in com- 
mon, both in the flare SOL2011-02-15 and in other flares, acoustically active or 
not: They all show strong, apparently transient signatures in vast regions from 
which no significant seismic emission emanates. It appears that there has to be 
something special about seismic-source regions such that any of the mechanisms 
yet considered (or perhaps others) could have a role in transient seismic emis- 
sion. For either chromospheric heating or Lorentz-force transients to satisfy the 
helioseismic signatures, something distinctive about their respective operations 
begs to be identified that makes it horizontally selective of the highly restrictive 
regions from which seismic radiation emanates. This is the primary theme of 
this study. Failing such an identification, it should begin to appear that seismic 
transient emission must depend upon some mechanism not yet considered. 

In the case of Lorentz-force transients, for example, a basis for "regional 
selectivity" might be the difHculty of credible magnetic measurements during the 
impulsive phases of flares. If transient magnetic signatures during the impulsive 
phase of the flare arc simply spurious, then magnetic variations that actually 
occur in the impulsive phase may be insufficiently sudden to excite a seismic 
transient — except where the signatures of transient seismic sources actually ap- 
pear. Magnetography is known to be highly susceptible to this liability under 
flaring conditions. If this is the basis of the regional selectivity we find, then 



^The spatial resolution of the HXR observations, at ~1.7 Mm. 
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50 Mm 



Figure 1. Doppler egression-power maps from SDO/HMI observations representing 1-mHz 
bandpasses centered at 3 mHz (upper left) and 6 mHz (upper riglit) in units of m^ s~'^. 
Middle row shows a pre-flare HMI continuum intensity snapshot (left) and an AIA snapshot 
(right) of 160 nm continuum radiation in the early impulsive phase, both normalized to unity 
for the quiet Sun. Bottom row shows a pre-flare HMI line-of-sight magnetogram (left), and 
line-of-sight difference in the magnetic field from before to after the impulsive phase, the time 
difference being 15 min, both in units of kGauss. Arrows locate compact sources of seismic 
emission in the upper right egression-power map. Concentric circles in the upper left frame 
show the dimensions of the pupil applied in the hclioseismic egression computations in the top 
row. 
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what is distinctive about seismic source regions so tliat magnetic transients in 
them are more sudden than elsewhere? 

For transients driven by impulsive heating, a highly suggestive basis for re- 
gional selectivity is the impact of a strong magnetic field on the overbearing 
radiative losses discussed in §1.1, which discourage the downward conduction 
of seismic transients in non-magnetic simulations. In the absence of a magnetic 
field, a shock driven by the expansion of the heated medium compresses the 
underlying medium, heating it and increasing its density. Both the heating 
and the increased density make the gas a more effective radiator. Hence, the 
radiative losses increase strongly, nonlinearly, with the degree of compression, 
which is 1-2 orders of magnitude in the low chromosphere in the simulations of 
Fisher, Canfield & McClymont (1985a, b,c). The onus appears to be whether any 
amount of impulsive heating can manifest a seismic transient that can survive 
the radiative losses entailed in passage through both the low chromosphere and 
photosphere to penetrate into the solar interior. 

As a screening measure to judge whether impulsive heating merits serious 
consideration, simulations such as those of lFisher et al.l(|2012p and lAUred et al.l 
(|2005p can seem discouraging indeed. However, we want to point out just the 
opposite, that these simulations are entirely consistent with helioseismic sig- 
natures as we know them. Impulsive chromospheric heating simulations to date 
have been computed only for non-magnetic media, from which detectable transient 
seismic emission has indeed never emanated, to our knowledge. Based upon this, 
we propose, as a hypothesis, that the very radiative losses that discourage the 
impulsive heating as a driver of transient seismic emission in the context of a 
non-magnetic medium are the key to the high degree of regional selectivity shown 
by transient seismic emission in flares. Essentially all transient seismic sources 
appear in sunspot penumbrae, regions known to be infused with strong, highly 
inclined magnetic fields. We want to propose that a major mechanical role of 
this field can be a radical reduction in radiative losses in those localities. This 
reduction in radiative losses could rescue impulsive heating as a tenable source of 
the seismic emission that emanates from some flares. Moreover, if backed up by 
appropriate MHD simulations that incorporate oblique magnetic fields, and by 
appropriate observational diagnostics, it may very well explain the high regional 
selectivity that characterizes transient seismic sources. 



3. The Mechanical Role of an Inclined Magnetic Field in 
Seismic- Transient Conduction 

For a rough, preliminary understanding of the role of an inclined field, we con- 
sider the simplest case, that of a horizontally uniform medium infused with a 
horizontal field, likewise horizontally uniform, into which we introduce a hor- 
izontally uniform shock, driven by a horizontally uniform thick-target-profile 
heating. In the MHD context the action of the horizontal magnetic field is 
essentially equivalent to that of a large increase in the compressional modulus 
of the medium, the magnetic component of which is 

KM = - = 2pA/, (1) 
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where 



B2 



(2) 




is the magnetic pressure, with B the magnetic field strength, and V represents 
the volume of the local unit of medium being vertically compressed by the shock. 
The compression, then, is perpendicular to B, the local magnetic induction. 
For a magnetic field of 500 Gausf|f| in the primary source region of the fiare 
SOLA2011-02-15, this is -^40 times the isothermal gas modulus. 



where Pg is the local gas pressure. 

The presence of a horizontal magnetic field whose strength is typical of sunspot| 
penumbrae, then, will make the medium much more "vertically rigid," i.e., 
much more resistant to compression perpendicular to the field. The medium 
will therefore conduct seismic energy downward with a much smaller degree of 
compression. The radiative losses will be reduced accordingly. Since the radiative 
losses in the non- magnetic medium increase strongly, i.e., highly non- linearly, 
with compression, wc have to suspect that the radiative losses will be a small 
fraction of those indicated by the non-magnetic simulations, for a sufficiently 
inclined field. This will introduce a high degree of selectivity towards regions of 
highly inclined magnetic field. In the ideal case of a horizontal field conducting 
a downwardly propagating transient, nearly all of the compressional work will 
be done on the magnetic field, which is almost completely elastic, hence lossless. 

The role of inclined magnetic fields, then, needs to be seriously considered as a 
prospective explanation of strong regional selectivity of transient seismic sources 
in favor of strongly inclined magnetic fields. Indeed, among the most conspicuous 
such source regions is the penumbral neutral line separating umbrae of oppos- 
ing polarities in (5-configuration sunspots. In these instances, e.g.. IMoradi et aP 
(|2007p (see Figure 2 of this study), the magnetic field is exceptionally strong, 
leading to an especially high magnetic modulus, fc^/, and is essentially horizon- 
tal, hence avoiding losses due to coupling between fast and slow modes (see 
ICally and HaSseSl [mT|) . 

4. Qualifications 

4.1. Magnetic Encumbrance of the Work Capacity of a Heated Gas 

Somewhat anticipating realistic simulations, it should be noted that the mag- 
netic rigidity that should (we propose) so greatly reduce radiative losses once a 



^This exercise selects the lower limit of our understanding of the typical range of the horizontal 
components of penumbral fields as an example. |Borrcro and Icl iimoto, 2011[ l find penumbral 
fields with horizontal components in the range 1000-1200 Gauss, with a stochastic scatter of 
±500 Gauss. 




(3) 
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Figure 2. Doppler egression-power map (lower right) of seismic emission from NOAA 
AR10720 at 00:42 UT, during the impulsive phase of the flare of 2005 January 15. The 
line-of-sight magnetogram (lower left) identifies AR10720 as a (5-configuration sunspot. Up- 
per left frame shows the pre-flare continuum intensity, from SOHO/MDl. Upper right frame 
shows the continuum-intensity excess in the impulsive phase of the flare, from GONG. Bottom 
left frame shows a line-of-sight magnetogram shortly before flare onset. Concentric circles in 
the upper left frame show the dimensions of the pupil applied in the helioseismic egression 
computation. 

transient has development must also inhibit this development. This is because it 
encumbers the expansion of the heated medium, which is the agent of mechanical 
work upon its surroundings. Heating applied to a low-/3, i.e. magnetically rigid, 
medium will respond with less expansion than a non-magnetic medium. At some 
point this will limit the capacity of the heated gas to do work on its surroundings, 
including on the magnetic field itself, hence the efficiency with which impulsive 
heating will drive a seismic transient. For an estimate of this encumbrance over 
a range of temperatures to which the gas can be heated, we determine the 
"mechanical work capacity" with and without a magnetic field for a unit mass 
(1 gram) of hydrogen, with a helium abundance of 10%, initially at temperature 
To of 4,000 K and under a gas pressure, Pgo, of 500 dyne cm~^. The gas is first 
heated at constant volume, Vq, to a temperature Tq -I- AT, raising its pressure 
as illustrated by the vertical displacements in the p-V trajectories plotted in 
Figure 3. In this representation, the non-magnetic gas is represented by the blue 
locus. The green locus represents an identical gas with a horizontal magnetic field 
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p— V Diagram 




10^ ]0'° 
V(cmVg) 



Figure 3. Trajectories in p-V space of a unit mass (1 g) of chromospheric hydrogen-helium 
initially at a temperature T of 4,000 K and pressure of 500 dyne cm"'' heated to 10^ K, then 
expanded adiabatically to the initial prc-hcated pressure, and finally cooled at that pressure. 
The blue locus represents a non-magnetic medium. The green locus represents the same cycle 
for a medium infused with a 500-Gauss magnetic field. We recognize the respective "adiabatic 
work capacities," AW, of these cycle to be the area in p-V-space enclosed by the respective 
curves. For the foregoing parameters, AW is f.44xf0^"' erg g"^ for the non-magnetic medium 
and 0.459x10^^ erg g~^ the magnetic medium. 

initially of 500 Gauss. Once heated, these gas samples are expanded vertically 
and adiabatically to their original pressures, and subsequently allowed to cool 
at the original pressure to the original temperature of 4,000 K. We recognize the 
areas of the respective closed curves in the p-V plane as the respective adiabatic 
work capacities, 

AW = / pdV, (4) 

of the heated gas samples. 
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Figure 4. Mechanical work capacity of an unmagnetizcd heated gas (blue) is compared with 
the same for a gas infused with a horizontal magnetic field whose initial strength is 500 Gauss 
is plotted as a function of the temperature increase, AT = Ti — Tq, in the heated gas. The 
gas, in thermal equilibrium with a helium abundance of 10%, is heated at constant volume 
from an initial temperature of 4,000 K to the temperature indicated by the abscissa. It is then 
expanded adiabatically until the composite of its gas and magnetic pressures has returned to 
the pre-heated composite. Ordinate values indicate the work done by the expanding gas and 
magnetic flux to occupy a greater volume than before the gas was heated. 

Figure 4 shows comparative plots of AW for temperature increments up 
to 10^ K. The capacity of a magnetized layer to do mechanical work on its 
surroundings is seen to roughly approach that of the non-magnetized gas when 
the heating is sufficient to increase the /? of the gas to something like unity or 
more. I.e., nowithstanding that the magnetic modulus, kmi strongly dominated 
gas dynamics in the low chromosphere in pre-flaring conditions, if the heating is 
sufficient to elevate the gas modulus, Kg, so that it dominates over the magnetic 
modulus, then the impact of the magnetic field upon the work capacity of the 
gas becomes relatively modest. If heated to coronal temperatures, i.e., 10^ K, a 
500-Gauss magnetic field reduces the work capacity by about a third. 

A magnetic medium, then, is considerably more conducive to seismic-transient 
generation if the medium heated is relatively dense. This suggests that higher- 
energy particles — which, in a thick-target model, for example, can penetrate 
to greater column densities — might be more conducive to seismic transients 
than low-energy particles, which would expend their heat in low-density media. 
This might explain the strong association of seismic transients with hard X-ray 
sources. 
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4.2. Regional Selectivity of Lorentz-Force Transients 

It should also be noted that seismic emission due to Lorentz-force transients, 
could be fairly selective in favor of regions of inclined magnetic field — under 
conditions that are similarly subject to observational control. This could be 
the result of magnetic measurements sampling a layer not very far above that of 
"unity beneath which the gas pressure, Pg, dominates the magnetic pressure, 
Pm, and, due to the high conductivity of the solar plasma, the magnetic field 
is rigidly frozen into it. HMI magnetic diagnostics, for example, are based upon 
the line Fe I 6170 A, whose core forms about 100 km above the height at which 
t(5000 A) passes through unity. The gas pressure, Pg, at this height is about a 
third of the magnetic pressure, pm, for a 500-Gauss field. Proceeding downward, 
from 140 km beneath this level, the gas pressure dominates, e-folding about 
every 100 km|f| The horizontal scale of seismic sources such as that labeled "SI" 
in Figure [1] is ^2.5 Mm, about 20 times the e-folding depth in pressure. 

We will now consider the simple model in which the magnetic flux we ob- 
serve passes from the surface in which it is measured into a rigidly frozen-in 
condition in an inflnitesimally thin layer beneath it, calling this the "thin- 
unity-/3-layer approximation," or simply the "thin-layer approximation" for short 
(see lAlvarado Gomez et a"Lll2012p . For simplicity, we represent Lorentz forces in 
the plane-parallel approximation adopted by [Hudson, Fisher and Welsch| (POOSP . 
The local transient pressure exerted upon the photosphere by a magnetic tran- 
sient, ^B, localized thereat can be expressed by 

SpM = ^{2B,6B, ~ B^-JB,,). (5) 

In the thin-layer approximation, the surface flux density, i.e., B^, is held constant 
during the flare, since the flux is rigidly frozen into the dense, highly conducting 
medium just beneath where it is observed. Hence, 5Bz in equation ([S]) must be 
null, whence 

= -^Bh-5Bh. (4) 

Under the thin-layer approximation, then, seismic transient excitation due 
to differential Lorentz-force transients, JB^, would vanish when B is vertical, 
lending a distinct preference to regions whose magnetic fields have a strong 
horizontal component. 

Whether the thin-layer approximation is a realistic representation of MHD 
in the environments of seismic sources to a tenable precision on the impulsive- 
phase time scales required for efficient transient seismic emission has not yet 
begun to be realistically addressed to our knowledge. The point of the foregoing 
exercise, then, is not to promote the accuracy of the thin-layer appoximation. 



*The foregoing parameters actually describe the "standard (quict-Sun) model" of 
[Christensen-Dalsgaard ProfEtt and Thompson| II1993II . For the sunspot penumbra, the tran- 
sition will have to be sharper, because of the Wilson-depressed photosphere of the sunspot 
l |Lindsey, Cally and Rempel, 2010[ l. 
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nor to encumber Lorentz-force-driven seismic transients with any of a range of 
qualifications that would follow from it. It is only to render a recognition of 
how some degree of regional selectivity could apply to seismic sources driven 
by Lorentz-force transients as we propose they do to impulsive chromospheric 
heating. The more important point is that we have abundant resources with 
which to address this question in lieu of having the answer to it now. 

5. What Should We Do? 

Considerable attention is now being paid to Lorentz-force transients as a prospec- 
tive source of transient seismic emission, motivated by |Hudson, Fisher and Welschj 
(|2008p and lFisher et al.| (j2012|) . This is only appropriate. If Lorentz-force tran- 
sients contribute significantly to transient seismic emission, the implications 
respecting the flare mechanism at large have to be formidable. If not, the im- 
plications of "not" arc probably similarly formidable. This might suggest an 
onset mechanism that, at least on the time scale of seismic transients, delivers 
a great deal of energy to a very small mass, i.e., high energy electrons, but, 
as IZarro and Canfieldl (|1989|) suggest, need not impart anything like the 10^^- 
dyne-s impulse required by a seismic wave. The primary goal of this study is to 
establish the need for a commensurate expenditure of attention on the contribu- 
tion of impulsive chromospheric heating, particularly the need to include inclined 
magnetic fields in simulations thereof. An appreciation of this distinction, if it 
applies, could lend us some very cogent clues that we have yet to secure about 
how flares work. 



6. How Do We Do It? 

Even 1-D simulations in a medium with an obliquely inclined fleld will require a 
considerable extension of standard non-magnetic 1-D simulations of impulsive- 
heating dynamics. The Lorentz force will likewise be oblique, introducing a 
horizontal component of motion. This will introduce the complication of cou- 
pling between modes whose restoring force is primarily compressional and modes 
whose restoring force is primarily magnetic tension ( jCally and Hansen] I201ip , 
the latter of which will add to the depletion of the former, hence of the eventual 
helioseismic signature — thereby further contributing to the regional selectivity 
attendant to inclined fields. The horizontal motion will break azimuthal symme- 
try, complicating the RT computations. In particular, the relative Doppler-shifts 
in the chromospheric lines will depend upon the azimuth of the direction of the 
radiation. So, line intensities, opacities, source functions, etc., will require an 
account for azimuthal dependencies. Realistic 3-D simulations will entail further 
complications still, and a greater expense in computing resources. However, for 
the special case of 1-D simulations of a downwardly propagating shock in a 
horizontally magnetized medium, i.e., a perpendicular wave or shock, azimuthal 
symmetry is restored, and the foregoing complications essentially vanish. This 
problem can be addressed by incorporating the magnetic modulus expressed by 
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equation (1) into the the gas dynamics as if this is simply a thermodynamic 
property of the gas. 

A feasible approach to the issue of selectivity, then, would be adaptations of 
existing HD algorithms applied to simulations such as those applied bv lFisher et aH 
(|2012p and lAUred et al.l (j2005p to a non-magnetic medium to represent a hori- 
zontally magnetized medium by the inclusion of an appropriate elastic modulus 
(equation 1) in addition to the thermally induced forces. If these computations 
indicate the survival of transient seismic emission penetrating into the solar 
interior consistent with hclioseismic signatures from sunspot penumbrae, this can 
serve as a practical basis for the investment needed for simulations of chromo- 
spheres in which the magnetic inclination is other than horizontal, i.e., oblique. 
These simulations will be needed for an assessment of the degree of regional 
selectivity that realistically characterizes impulsive chromosphcric heating as a 
contributor to transient seismic emission as indicated by helioseismic signatures. 
At the very least, it is important to recognize that the relatively simple extension 
of the 1-D simulations of IFisher et all (PUT^ lAllred et all (^005)) to include 
horizontal magnetic fields is very timely at this juncture. Given the strongly 
selective character of transient seismic emission in regions of strongly inclined 
magnetic field, impulsive chromosphcric heating as a driver of transient seismic 
emission should not be ruled out until at least this adaptation has been made 
and applied at length. 

7. Summary 

Simulations of impulsive chromospheric heating as a driver of transient seismic 
emission from flares show strong chromospheric shocks with energies roughly 
consistent with helioseismic observations of this phenomenon. However, these 
simulations also show radiative losses that deplete nearly all of the energy 
initially investmented into the shock as it plows downward into the underlying 
chromosphere and photosphere. This has introduced doubt that sufficient energy 
from such a shock can penetrate into the solar interior to match that implied 
by the helioseismic signatures. The inclusion of a highly inclined magnetic field, 
typical of transient-seismic-source regions would have a major impact upon this 
assessment. The magnetic field can greatly increase the compressional modulus 
of the chromospheric medium, greatly reducing the compression, hence the ra- 
diative losses that act once a shock has developed. This might explain the strong 
regional selectivity transient seismic sources show favoring sunspot penumbrae, 
where magnetic fields are highly inclined. This includes the frequent occurrence 
of seismic emission from the magnetic neutral lines separating opposing polarities 
in ^-configuration sunspots. 

On the other hand, the infusion of a plasma with a strong magnetic field can 
inhibit the work capacity of such a medium in response to impulsive heating. 
Because of this, transient generation due to impulsive heating of a magnetic 
medium favors denser media. This suggests further conditions for regional selec- 
tivity, which might relate to the affinity of transient seismic sources to regions 
from which hard X-rays are seen to emanate. 
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In the context of our present knowledge, the foregoing hypotheses are highly 
speculative. However, the crucial point is that the inclusion of inclined magnetic 
fields is fundamental to our understanding of the role of impulsive atmospheric 
heating in transient seismic emission. 

The introduction of obliquely inclined fields will complicate the HD simu- 
lations considerably. However, horizontal magnetic fields can be incorporated 
into existing in 1-D HD codes relatively simply. The results of relatively cheap 
simulations such these can serve as a practical basis for more realistic, but 
expensive, 1-3-D simulations that incorporate oblique or stochastic magnetic 
fields. 

Bear in mind that the purpose of this study is not to promote atmospheric 
heating of any kind as "the source" of transient seismic emisson in acoustically 
active fiarcs. It is rather to establish the essentiality of a realistic account for 
inclined magnetic fields in simulations of this phenomenon for a realistic assess- 
ment of its contribution to transient seismic emission. Nor should we relax into 
the supposition that impulsive atmospheric heating is by any stretch the only 
agent that offers the regional selectivity shown by the helioseismic signatures. 
It simply must be understood that inclined magnetic fields arc fundamental to 
the contribution of impulsive atmospheric heating as a realistically prospective 
contributor to transient seismic emission from fiarcs. This mechanism should 
not be ruled out of consideration as a major contributor, even possibly to 
sole contributor, to transient seismic emission without this adaptation of HD 
simulations of this phenomenon having been developed and applied at length. 
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